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AN  AUTOMATED  MEASURING  SYSTEM  FOR  ELECTROACOUSTIC  TRANSDUCERS 


INTRODUCTION 

The  precise  measurement  of  immittamce  is  a  valuable  tool  for  the 
analysis  of  electroacoustic  transducers  and  piezoelectric  materials. 
Imnittance  is  used  in  the  determination  of  transducer  efficiency, 
electrical  matching,  and  the  tuning  and  analysis  of  transducer  performance 
and  of  material  parameters  of  active  transducer  materials.  Previous 
techniques  for  immittance  measurements  involved  the  use  of  such  devices  as 
impedance  bridges,  analog  voltmeters,  phase  meters,  and  some  special  analog 
devices  like  the  vector  impedance  meter.  All  of  the  before  mentioned 
devices  suffer  from  one  or  more  restrictions;  e.g.,  speed,  precision,  or 
stability. 

The  HP  Model  4192A  Digital  Network  Analyzer  is  computer  controlled 
and  has  the  ability  to  make  immittance  measurements  at  high  speed  with 
excellent  precision  and  stability.  Combined  with  the  developed  circular 
approximation  algorithm,  the  measurements  cam  be  made  with  greater 
efficiency  to  provide  the  maximum  amount  of  information  with  the  minimum 
number  of  data  points.  This  method  also  allows  for  fast,  accurate 
computation  of  critical  points  of  immittamce  without  knowing  the  series 
(parallel)  resonamce  of  the  transducer.  The  optimization  of  data  points  is 
achieved  by  the  determination  of  frequencies  that  will  be  quasi-uniformly 
distributed  axound  the  immittamce  circle  as  opposed  to  the  "clustering"  of 
frequencies  above  and  below  the  frequency  of  maximum  immittamce. 

Figure  1  compares  am  immittamce  circle  defined  by  a  lineair  sequence  of 
frequencies  with  one  in  which  the  frequencies  are  generated  from  the 
circular  approximation  method.  There  axe  two  hundred  points  in  the  lineax 
frequency  sweep  of  Fig.  la  and  only  fifty  points  in  Fig.  lb.  This  circular 
approximation  method  has  been  implemented  on  various  piezoelectric 
tramsducers  amd  works  exceptionally  well  in  most  cases.  The  timesaving  is 
greatest  when  measuring  high-Q  tramsducers.  Linear-spaced  frequency  points 
appear  diluted  on  the  steep  slopes  of  the  imaginary  part  of  a  high-Q 
immittamce,  as  shown  in  Fig.  2.  For  low-Q  curves  the  frequency 
distribution  is  more  evenly  spaced.  When  the  data  are  presented  in  the 
form  of  Imaginary  verses  real  immittamce,  the  data  from  high-Q  systems  are 
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Fig. la  -  Immittance  loop  showing 
the  distribution  of  data 
points  resulting  from  a 
linear  sweep  of  two 
hundred  frequency  steps. 


Fig.  1b  -  Immittance  loop  showing  the  more 
uniform  distribution  of  only  50 
data  points  when  frequencies 
are  determined  by  the  circular 
approximation  method. 


Fig.  2  -  Distribution  of  data  points  for  equally  spaced  frequencies  on  the  imaginary 
part  of  the  immittance  curve  for  a  high-Q  system  and  for  a  low-Q  system. 
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concentrated  near  the  beginning  and  end  of  the  iinmittance  circle.  For  low- 
Q  systems  the  data  points  approach  a  uniform  distribution. 

CIRCULAR  APPROXIMATION  ALGORITHM 

The  circulax  approximation  method  can  be  applied  to  admittance  and 
impedance  loops.  Because  of  the  similarity  of  application  to  impedance  and 
admittance,  and  to  avoid  writing  "impedance  or  admittance",  the  two 
quantities  axe  implied  when  the  term  Immittance  is  used  in  the  following 
text.  The  differences  between  the  admittance  terms  and  the  impedance  terms 
will  be  discussed  in  the  PROGRAM  section  of  this  report.  Thoughout  this 

report  the  word  "points"  appears  many  times  with  vaxious  meanings.  To 

assure  clarity  three  symbols  are  defined: 

o  The  term  "points ( f)"  refers  to  a  frequency. 

o  The  term  "points(c)"  refers  to  the  ordered  pair 

on  the  the  theoretical  circle. 

o  The  term  "points(i)"  refers  to  the  real  and 
imaginary  position  on  the  immittance  loop. 

* 

The  method  presented  here  will  be  based  on  the  following  equation: 

lL2  -  U.2] 

^  — ■  C1] 

This  equation  is  derived  [1]  from  a  simple  equivalent  circuit  of  a 
piezoelectric  resonator,  where  L  is  the  inductance,  o>  is  the  angular 
frequency,  c«)i  is  a  resonant  frequency  (series  resonance  when  measuring 
admittance,  parallel  resonance  when  measuring  impedance),  ©  is  the  phase 
angle,  and  R  is  the  resistance.  Cady  uses  u)Q  (angular  frequency  of  zero 
crossing)  in  Eq.  (1)  instead  of  u>^.  But  is  also  an  angular  frequency  at 
a  zero  crossing  since  a  transformation  of  axes  is  performed  with  the 
circulax  approximation  method.  It  is  algebraically  convenient  to  write 
Eq.  (1)  in  the  following  form:  Given  u>,  R,  and  ©  for  three  points(f), 
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and  L  are  determined  by  applying  Cramer's  rule  [2]  to  this  expression i 

2  1  2 

a>i  +  —  o)Rtan(G)  «  u>  .  [2] 

By  using  the  following  equation 

2  2 

Lu>  -  Lok  -  G>Rtan(e)  «  0,  [3] 

the  angular  frequencies  (w)  are  determined  as  a  function  of  the  angle  ©  by 
solving  the  quadratic  equation.  Equations  (2)  and  (3)  axe  the  key 
equations  used  in  the  circle  approximation  method  and  will  be  referred  to 
many  times  in  the  text  of  this  paper.  The  argument  used  for  applying 
Eqs .  ( 1 )  through  ( 3 )  to  both  impedance  and  admittance  is  found  in  Appendix 
A. 


The  algorithm  used  in  determining  the  admittance  or  impedance  loops 
consists  of  two  major  parts.  First,  the  "best  fit"  circle  that 
approximates  the  inmittance  loop  is  determined  from  parameters  of  the 
simple  equivalent  circuit  of  a  piezoelectric  resonator.  Secondly,  the 
circular  approximation  is  then  divided  into  equal  arcs  that  are  used  to 
determine  quas i-equally  spaced  points(  i )  on  the  immittance  loop. 

To  determine  the  "best-fit"  circle,  one  must  first  identify  any  three 
frequency  points(f)  on  the  immittance  loop  with  the  abscissa  and  ordinate 
identified  as  conductance  and  susceptance  (or  resistance  and  reactance), 
respectively.  If  previous  knowledge  exists  about  the  transducer  to  be 
measured,  these  three  points(f)  are  easy  to  identify.  If  however,  there  is 
no  previous  knowledge  about  the  transducer,  the  program  will  assist  the 
user  in  determining  the  three  points(f). 

The  HP  4192  Impedance  Analyzer  measures  the  real  and  imaginary  parts 
of  initiittance  at  each  of  these  frequency  points(f)  and  the  ordered  pairs  of 
data  are  stored  as  (x^,y^),  where  i  goes  from  1  to  3  corresponding  to  each 
of  the  three  frequency  points(f).  These  ordered  pairs  are  used  in  the 
general  form  of  the  equation  of  a  circles 
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X±D  +  y±E  +  F  -  -[x±2  +  y.2] 


C4] 


of  center 


and  radius 


+  E 


,1/2 


4Fj 


C5] 


To  find  D,  E,  and  F,  we  impose  the  condition  that  the  coordinates  of  each 
of  the  given  points(ific)  satisfy  Eq.  (4).  When  these  values  for  D,  E,  and 
F  are  substituted  in  Eq.  (5),  we  have  the  radius  and  center  of  the  circle 
through  the  given  points(ific).  By  using  the  calculated  radius  and  center, 
two  points(  c )  located  120  degrees  from  a  fixed  initial  point( i&c )  are  found 
as  shown  by  the  boxes  in  Fig.  1.  The  authors  have  arbitrarily  chosen  the 
fixed  point(c)  to  be  the  second  initial  point(  i ) .  One  may  use  any  initial 
point(i)  as  long  as  it  remains  consistent  throughout  the  entire  procedure. 


In  order  that  this  procedure  apply  both  to  loaded  and  unloaded 
transducers,  a  translation  of  axes  is  performed  on  all  points( i )  used  in 
Eqs .  (2)  and  (3).  The  formulas  for  the  translation  of  axes  are: 


x  =  x  -  x  +  r  and 

to  c 


yo  - 


[6] 


where  (xc  +  r,yc)  is  the  new  origin.  The  ordered  pair  (xc/Yc)  is  the 
center  of  the  circle,  r  is  the  radius,  (x0/Y0)  are  the  old  coordinates,  and 
(xt,yt)  are  the  old  coordinates  with  the  new  origin.  The  slope  [tan(0)]  is 
determined  by  first  applying  Eqs.  (6)  to  a  given  point(i).  The  resistance 
R,  also  found  in  the  Eqs.  (2)  and  (3),  is  replaced  by  the  diameter  or 
reciprocal  of  the  diameter  of  the  circle  depending  on  which  immittance  loop 
is  chosen.  The  quantities  and  L  are  calculated  using  Eq.  (2)  and  the 
data  of  two  intial  points(  i )  along  with  information  from  the  circle.  With 
the  calculated  parameters  (to^,  L,  and  ©)  of  the  two  120-degree  points(c), 
the  frequencies  (o)/2tt)  at  the  120-degree  points(c)  are  calculated  by 
applying  Eq.  (3).  The  HP  4192  Impedance  Analyzer  measures  the  conductance 
and  susceptance  (or  resistance  and  react ance )  at  each  frequency,  and  the 
points(  i )  are  stored.  These  points( i )  are  illustrated  by  the  "triangles" 
in  Fig.  3 . 
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Fig.  3  -  An  immittance  loop  with  a  first-circle  approximation. 

Dots  indicate  the  first  three  frequencies  chosen. 

Squares  mark  the  120°  points. 

Triangles  represent  associated  points  on  immittance  loop. 


Now  a  new  circle  is  calculated  using  the  two  120-degree  points(c) 

(triangles  in  Fig.  3)  and  the  second  intial  point(i&c).  Note  that  this  new 
circle  will  be  a  better  approximation  to  the  immittance  loop  than  was  the 
previous  one.  This  procedure  iterates  until  two  successive  circles  are 
within  the  desired  tolerance.  The  result  of  this  convergence  is  a  "best- 
fit”  circle,  which  is  then  divided  into  equal  sectors,  as  shown  in  Fig.  4; 
the  frequencies  are  found  at  each  "triangle”  location  by  applying  Eq.  (3). 


Fig.  4  -  An  immittance  loop  with  a  first-circle  approximation. 

Equally  spaced  points  on  the  circle  create  quasi-equally 
spaced  points  on  the  immittance  loop. 

The  HP  4192  Impedance  Analyzer  measures  the  conductance  and  susceptance  at 
the  "triangle”  frequencies,  and  the  immittance  values  are  stored  in  a 
file.  Note  the  dots  on  the  immittance  loop  in  Fig.  4  are  not  exactly 
"equally  spaced",  but  certainly  there  is  a  better  distribution  of  pomts(i) 
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than  is  obtained  from  a  linear  frequency  sweep.  This  is  what  we  previously 
referred  to  as  quasi-equally  spaced  points(i). 

CRITICAL  FREQUENCIES 

Once  the  loop  is  characterized  with  the  desired  number  of  quasi- 
equally  spaced  points(i),  specific  frequency  points(f)  are  needed  for  later 
calculations  of  various  transducer  parameters.  The  frequency  points(f)  [3] 
that  are  of  interest  are: 


V 

frequency  at  maximum  admittance  (minimum  impedance) 

fn‘ 

frequency  at  maximum  impedance  (minimum  admittance) 

V 

frequency  at  maximum  conductance 

V 

frequency  at  maximum  resistance 

fr: 

resonance  frequency  ( susceptance  is  zero ) 

fa‘ 

antiresonance  frequency  (susceptance  is  zero) 

fh  ! 
a 

frequency  at  maximum  susceptance 

fi  : 
a 

frequency  at  minimum  susceptance 

fh  •• 
i 

frequency  at  maximum  reactance 

fi  ! 
i 

frequency  at  minimum  react ance 

When  characterizing  an  admittance  loop,  the  frequency  points(f)  fm, 
fg,  fr'  fp'  ^haf  fla  axe  determined  by  a  maze  technique  that  will  be 

described  in  the  program  section  of  this  report.  When  an  impedance  loop  is 
characterized,  the  frequency  points(f)  fa,  fp,  fn,  fs,  f^,  and  f ^  are 
determined.  One  may  note  that  the  frequency  sets  differ  when  calculating 
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impedance  or  admittance  loops.  The  authors  believe  that  the  two  most 

important  points(f)  are  the  series  and  parallel  resonance  ( f g  and  fp). 

These  points(f)  can  easily  be  found,  and  they  occur  for  both  loaded  and 

unloaded  transducers.  Both  of  these  conditions  do  not  apply  to  the 

frequency  pairs  fm,  fn  and  fa/  fr*  To  determine  the  frequencies  fm,  fn  one 

must  do  a  frequency  sweep,  which  can  be  rather  time  consuming.  The 

points(f)  fa/fr  are  easy  to  find  by  implementing  root-finding  techniques; 

however  these  points(f)  do  not  always  occur  in  loaded  transducers.  For 

these  reasons,  the  focus  has  been  on  obtaining  f0,f_.  The  reader  may  note 

s  F 

that  fQ  is  the  frequency  at  maximum  conductance;  however,  it  is  not  the 
frequency  of  minimum  resistance  as  can  be  seen  in  Fig.  5. 


UJ 

Z 

< 
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UJ 
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CONDUCTANCE 


Fig.  5  -  Impedance  loop  (5a)  and  admittance  loop  (5b), 
showing  their  critical  frequencies. 

Analogously,  fp  is  the  frequency  at  maximum  resistance;  but  it  is  not  the 
frequency  of  minimum  conductance. 

The  points(f)  of  secondary  importance  axe  the  frequencies  of  minimum 

and  maximum  susceptance  ( f y.  and  f-,  ).  These  points(f),  sometimes  called 

na  x  a 

half-power  points(f),  are  used  to  compute  the  electrical  quality  factor  Q. 
We  also  compute  the  frequencies  of  minimum  and  maximum  reactance  ( f ^  and 
fj^).  These  points(f),  however,  are  not  the  half-power  points(f)  and  cannot 
be  used  for  calculating  the  quality  factor  Q.  The  only  purpose  for 
determining  these  points(f)  is  that  they  locate  the  frequencies  at  the  top 
and  bottom  of  the  loop.  The  remaining  frequencies  axe  easily  determined 
for  the  particular  immittance  loop  calculated,  and  they  serve  as  helpful 
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parameters  for  detecting  problems. 


PROGRAM 


The  computer  program  MI CAM,  presented  herein,  is  by  no  means  the  best 
implementation  of  the  circular  approximation  method.  However,  the  authors 
feel  confident  that  this  program  can  be  refined  in  a  manner  that  would  be 
applicable  in  any  research  or  manufacturing  environment.  The  MICAM 
program,  written  in  FORTRAN  77,  is  run  on  a  Digital  Equipment  Corporation 
VAX  11/780  computer  at  the  Underwater  Sound  Reference  Detachment  of  the 
Naval  Research  Laboratory  as  a  research  tool  as  well  as  a  quality-control 
problem  detector.  A  version  of  this  program  is  also  written  in  HP  Basic 
and  can  be  executed  on  the  HP  9836.  The  MICAM  program  has  been  tested  on 
various  transducers,  both  in  air  and  in  water,  and  seems  to  perform 
exceptionally  well.  It  does  have  limitations,  which  will  be  explained  in 
detail  in  the  RESTRICTIONS  section  of  this  report.  Since  the  program  is 
quite  involved,  the  line-by-line  descriptions  of  MICAM  and  the  program 
listings  are  included  in  Appendix  B.  The  purpose  of  examining  individual 
program  lines  is  to: 

o  Show  how  the  algorithm  has  been  implemented. 

o  Give  differences  in  admittance  and  impedance  loops. 

o  Explain  how  critical  frequencies  are  determined. 

The  program  MICAM  calls  three  subroutines;  SIMEQ,  RSORT,  and  FILE.  The 
SIMEQ  subroutine  solves  simultaneous  equations.  The  RSORT  subroutine  sorts 
numbers  in  descending  or  ascending  order.  And  the  FILE  subroutine  stores 
the  immittance  information  in  FORTRAN  files. 

RESTRICTIONS 


Limitations  exist  in  all  algorithms  and  should  be  used  as  one 
criteria  for  choosing  one  method  over  another.  The  following  paragraphs 
identify  two  examples  where  a  linear  frequency  sweep  would  be  preferable  to 
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determining  frequencies  by  the  circular  approximation  method. 

One  problem  occurs  when  running  very  low-Q  or  highly  damped 
transducers  where  no  loops  axe  formed.  Since  the  circular  approximation 
method  is  designed  to  emulate  the  inmittance  loop  with  a  circle,  the 
algorithm  fails  to  converge  on  very  low -Q  or  highly  damped  transducers. 

The  MI CAM  program  accommodates  this  problem  by  allowing  the  user  to  perform 
a  frequency  sweep.  If  the  transducer  is  of  low  Q,  a  linear  frequency  sweep 
will  give  quasi-equally  spaced  points(i),  a s  illustrated  in  Pig.  2. 

A  second  problem  occurs  when  two  or  more  resonances  are  in  the  same 
frequency  region.  They  may  appear  as  small  or  large  satellite  loops. 
Usually  the  small  satellite  loops  are  not  a  problem  and  MICAM  will  converge 
to  the  larger  loop.  Nevertheless,  when  immittance  is  the  combination  of 
loops  of  relatively  the  same  size,  MICAM  converges  to  one  of  them, 
depending  on  the  initial  frequency  given.  This  loop  may  not  necessarily  be 
the  one  desired.  However,  the  authors  feel  that  MICAM  may  be  optimized  for 
specific  applications  by  relatively  minor  programming  changes. 


CONCLUSION 


The  circular  approximation  algorithm,  when  used  with  the  HP  4192A 
Impedance  Analyzer,  makes  a  highly  efficient  system  for  measuring  impedance 
and  admittance  of  resonant  electromechanical  devices.  The  system  hats 
proven  to  be  applicable  to  the  measurement  of  active  electromechanical 
materials  as  well  as  composite  transducers.  The  algorithm  has  been  used 
with  excellent  results  on  a  VAX  computer  as  well  as  on  personal  computers. 

Although  specifically  developed  as  a  research  tool,  the  automated 
immittance  measuring  system  is  extremely  well  suited  for  quality-control 
and  production  testing  of  electromechanical  transducers  and  materials. 
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APPENDIX  A 


Analysis  for  Applying  Equations  to  Both  Inmittance  loops 

The  charact eristic  electrical  property  of  a  piezoelectric  resonator 
is  the  equivalent  series  chain  RLC  in  Fig.  Al.  The  graphical 
representation  of  the  RLC  branch  is  more  easily  represented  by  the 
admittance  which  is  represented  as  a  circle  referred  to  as  the  fiducial 
circle  [1].  The  graph  can  be  easily  extended  to  more  complicated  networks 
connected  in  parallel,  such  as  CQ. 


Fig.  Al  -  Simple  equivalent  circuit  of  a  piezoelectric  resonator. 


The  simple  form  of  admittance  of  the  RLC  chain  is  all  that  is  necessary  to 
approximate  the  fiducial  circle.  The  argument  extends  to  the  impedance 
representation  in  that  all  the  information  for  the  impedance  circle  is 
contained  in  the  admittance  circle.  If 

Y(u>)  =  G(u>)  -  jB(u>), 


then 


tan(0) 


G(  t*>  >  * 


[Al] 


Since  Z(w)  is  the  complex  reciprocal  of  Y(u>),  Z(w)  can  be  written  as 


Z(w) 


G(m)  +  B(o>) 

G2(w)  4-  B2(w)  G2(u>)  +  B2(w) 


[A2] 
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or  Z(  to )  cam  be  written  ate 


Z( (i>)  -  R(to)  +  jX( to) .  [A3] 

Equating  the  real  amd  imaginary  parts  of  Eqs.  (A2)  and  (A3)  results  in  the 
following  expressions: 


R(to) 


G(o>) 


2  2 
G  (to)  +  B  (to) 


and 


X(to) 


B(to) 

G2(  to  )  +  B2(g>) 


[A4] 


[A5] 


For  an  impedance  loop,  tan(0)  can  be  written  as 


tan(G) 


X(u) 

R(io)' 


[A6] 


or  can  be  equivalently  expressed  by  using  Eqs.  (A4)  and  (A5)  as 


tan(Q) 


-B(to) 
G(  to)* 


Note  that  Eqs.  (Al)  and  (A6)  axe  identical  except  for  being  opposite  in 
sign.  Therefore  for  any  to  in  the  admittance  circle  there  is  a 
corresponding  to  of  the  same  value  in  the  impedance  circle  with  the  same  but 
negative  phetse  value. 
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Program  Description  and  Listing 


DESCRIPTION 


This  appendix  contains  a  detailed  description  of  MI CAM  as  well  as 
program  listing.  The  program  description  supplements  the  CIRCULAR 
APPROXIMATION  ALGORITHM  section  of  this  report  and  is  necessary  for  a  full 
understanding  of  this  method.  Program  lines  1-76  are  array  declarations 
while  lines  94-145  contain  the  interactive  input  section.  There  are 
several  items  of  interest  in  this  section.  A  question  regarding  speed  is 
asked  at  line  98.  Low  speed  provides  an  average  measurement  mode 
(approximately  one  measurement  per  second)  to  obtain  measurement  data  of 
higher  resolution  and  repeatability  then  at  medium  speed  ( 5  measurements 
per  second)  or  high  speed  (10  measurements  per  second).  At  line  112,  the 
number  of  points (  i )  refers  to  the  number  of  quasi-equal ly  spaced  points( i ) 
that  axe  desired  on  the  immittance  loop.  In  lines  118  and  120,  the 
parallel  or  series  resonance  is  required  depending  upon  which  loop  is 
calculated.  These  frequencies  are  estimated  from  a  priori  knowledge.  The 
more  confident  one  is  of  these  points(f),  the  smaller  the  frequency  step 
size  can  be  taken  (line  125). 

Lines  138-163  is  a  section  of  code  that  allows  the  user  to  perform  a 
frequency  sweep.  This  is  executed  only  if,  for  some  reason,  the  circulax 
approximation  method  does  not  apply.  The  portion  of  code  that  includes 
lines  176-215  is  the  location  where  the  conductance  and  susceptance 
(resistance  and  reactance)  for  the  three  initial  points(f)  are  determined. 
The  ARRAY  variable  contains  the  coefficents  of  Eq.  4.  In  line  216,  the 
subroutine  SIMEQ  solves  the  three  simultaneous  equations,  which  ultimately 
results  in  the  center  and  radius  of  the  circle  given  in  lines  218-220. 

The  code  from  lines  221-237  is  used  to  solve  Eq.  (2)  for  and  1/L, 
which  are  used  to  calculate  the  frequencies  at  two  initial  points(i)  and 
the  motional  inductance  value  L  ( lines  241-250 ) .  Line  252  is  a  decision 
statement  that  determines  if  the  difference  between  two  successive  series 
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(parallel)  resonant  frequencies  is  within  a  desired  tolerance.  If  the 
frequency  is  within  tolerance,  the  program  has  determined  the  best  fit 
and  L  for  Eq.  (1).  If  the  tolerance  is  not  met,  a  new  and  hopefully 
better-fit  circle  will  be  computed  using  the  frequencies  at  the  two  120- 
degree  points(c)  (lines  261-283).  These  120-degree  frequency  points(f)  are 
found  by  solving  Eq.  (3)  and  substituting  the  recently  computed  values  for 
0)^  and  L.  This  entire  procedure  repeats  until  the  tolerance  condition  is 
met  in  line  252. 


The  second  major  computational  portion  of  the  MICAM  program  is  found 

in  lines  294—321.  This  portion  uses  the  and  L  that  meet  the  tolerance 

specifications  and  Eq.  ( 3 )  to  calculate  the  quasi-equally  spaced 

points(i).  The  remaining  portion  of  the  program  determines  the  critical 

frequencies  described  earlier  in  this  report.  The  RSORT  subroutine  sorts 

the  data  from  the  largest  to  the  smallest  value.  In  line  322,  the 

conductance  (resistance)  is  sorted  so  that  the  first  value  in  the  array 

will  be  the  f0  (f_).  The  RSORT  subroutine  is  called  again  in  line  327 
s  P 

where  frequency  is  the  variable  sorted.  The  remaining  critical  frequencies 
are  determined  by  exercizing  a  maze  technique.  This  technique  is  found  in 
lines  337-384  and  is  based  on  previous  knowledge  on  all  immittance  loops, 
shown  in  Fig.  5.  The  criteria  used  are: 


f 

m 


<  f 

8 


<  f 

r 


<  f 

a 


<  f 


<  f  . 
n 


[Bl] 


A  critical  frequency  point(  f )  is  searched  only  in  the  region  wherein  it 
could  occur;  i.e,  it  satisfies  the  inequalities  of  Eq.  (Bl).  All  of  the 
critical  frequencies  axe  chosen  from  the  points(ifif)  of  the  inmittance  loop 
except  for  fp(fs).  The  points(f)  fp(fs)  require  that  more  points(f)  be 
calculated  in  order  to  have  relatively  accurate  data.  For  example,  when  an 
admittance  loop  hats  been  run,  fp  is  determined  by  switching  the  HP  analyzer 
to  the  impedance  mode,  spotting  the  point(f)  f1  ,  and  sweeping  the 
frequency  until  the  frequency  at  maximum  impedance  is  found*  On  the  other 
hand,  when  an  impedance  loop  has  been  run,  fs  is  determined  by  switching 
the  analyzer  to  the  admittance  mode,  spotting  the  point(f)  ,  and  then 
sweeping  the  frequency  in  reverse  until  the  frequency  of  maxi  mum  admittance 
is  found.  This  technique  is  illustrated  in  lines  408-431. 
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The  remaining  portion  of  the  program  can  be  considered  a  formatted 
output  section  that  prints  a  table  of  all  critical  frequencies  described 
and  files  the  data  points(i)  for  plotting.  The  plots  that  can  be  created 
include  a  graph  of  magnitude  of  immittance  verses  frequency  and  a  graph  of 
the  imaginary  part  of  immittance  verses  the  real  part  of  immittance .  Also 
the  critical  frequencies  are  stored  so  they  may  be  overlayed  on  the 
immittance  plot  as  shown  in  Figs.  B1  and  B2,  which  are  examples  of 
admittance  loops  obtained  on  a  USRD  Type  F42A  transducer  in  air  and  in  the 
USRD  Lake  Facility  by  MI  CAM. 


Fig.  B1  -  Admittance  loop  obtained  on 
an  F42A  transducer  in  air. 


Fig.  B2  -  Admittance  ioop  obtained  on 
an  F42A  transducer  in 
the  USRD  Lake  Facility. 


There  are  several  differences  in  MICAM  for  impedance  or  admittance 
loops.  The  first  difference/  of  course,  is  that  the  HP  analyzer  is  in  the 
proper  mode  for  the  desired  measurements  (lines  171  and  172).  Secondly, 
there  is  a  difference  in  the  motional  resistance  (R)  found  in  Eqs.(l) 
through  (3)  for  the  two  loops.  In  the  admittance  mode,  R  is  approximated 
by  the  reciprocal  of  the  diameter,  and  in  the  impedance  mode,  R  is 
approximated  by  simply  using  the  diameter  of  the  circle.  Two  interactive 
runs  of  the  program  MICAM  are  included  in  Appendix  B,  which  illustrate  how 
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the  user  caui  obtain  results  in  either  iirmittance  node. 
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LISTING 


0001 

0002 

0003 

C 

MICAM  PROGRAM 

0004 

0005 

C 

WRITTEN  BY  TINA  RUGGIERO 

0006 

C 

FEBRUARY,  1983 

0007 

0008  C  THIS 
0009  C 


0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 


C 

C 


PROGRAM  DETERMINES  A  FREQUENCY  SWEEPING  FUNCTION  BY  APPROXIMATING 
AN  ADMITTANCE  LOOP  WITH  A  PERFECT  CIRCLE.  ALSO  FN,  FM, 

FAf  FR,  AND  ADMITTANCE  LOOPS  ARE  STORED  AND  ARE  AVAILABLE 
FOR  PLOTTING. 


DOUBLE  PRECISION  DEF  (3) ,  CONS (3) ,  ARRAY  (3, 3) ,  COEFF  (3, 3) ,  COLUMN  (3) ,  ANS  (3) 
REAL  RES(  1000) ,  XES(IOOO),  FREQ  ( 1000)  t  X  (6) , Y(6) ,  F(6)  ,  FR(3),  FRE(IOOO) 

REAL  K33,  CE,  CP,  COB,  CAP,  XCOF (5) ,  C0F(5) ,  ROOTR  (5) ,  ROOTI  (5) ,  ADM  (1000) 

REAL  CRX(8) ,  CRY (9) ,  MAG  (1000) 

COMMON  FR,  AMR,  RR,  XR,  FA,  AMA,  RA,  XA,  FM,  AMM,  RM,  XM,  FN,  AMN,  RN, XN,  XH,  XL 
COMMON  RH,  RL,  FH,  FL 
COMMON  /LENGTH/LEN 


BYTE  ANSWER  (80) 


BYTE  A(13) 

DATA  A  /'O',  'L',  0, 0, 0, 0, 0, 0, 0, 0, 0,  'E',  ’N7 
BYTE  B(12) 

DATA  B  /T, ’F\  0,0, 0,0,0,  O',  0,0,'E\'N7 
BYTE  C(13) 

DATA  C  /'PVF\0,0,0,0,0,0,0,0,0,'EVN7 
BYTE  D  (12) 

DATA  D  /'S’,  rF\  0, 0, 0,  0, 0, 0, 0,  0,  'E\  fN7 
BYTE  E (13) 

DATA  E  /'F\’R\0,0,0,0,0,0,0,0,0,'EVN7 
BYTE  IDSP  (4) 

DATA  IDSP  /’F',  rV,  'Af,  f47 

BYTE  IADV (2) 

DATA  IADV  /lW,,,27 

BYTE  IREV (2) 

DATA  IREV  /,W,,f47 

BYTE  IAD1  (2) 

DATA  I  ADI  f'U\'Vf 

BYTE  IAD0(2) 


I  VOLTS 

!  START 

I  STOP 

I  STEP 

I  SPOT 

D I  S 

{STEP  UP 

ISTEP  DOWN 


IDATA  READY  ON 


IDISPLAY  CAP 
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0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0060 

0081 

0062 

0063 

0064 

0085 

0066 

0087 

0068 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 
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DATA 

IAD0  /,D\’07 

DATA  READY  OFF 

BYTE 

IADM  (4) 

DATA 

lADM/'F',  'V,  'A*,  *2V 

DISPLAY  G 

BYTE 

IAV1  (2) 

(AVERAGE  ON 

DATA 

IAV1  /V,  *1V 

BYTE 

IAH1  (2) 

IHIGH  SPEED  ON 

DATA 

lAHI/'H'.'IV 

BYTE 

IAH0(2) 

IHIGH  SPEED  OFF 

DATA 

lAHO/'HVOV 

BYTE 

IAMM  (2) 

IMANUAL  SWEEP 

DATA 

lAMM/’W’/O'/ 

BYTE  IAC3(2) 

DATA  IAC3/*C\  '3'/ 

BYTE  IAC2C2) 

DATA  IAC2/'C\'27 

BYTE  IAEXC2) 

DATA  IAEX/*E\  *XV 

IADDR=1 
ITIMO=30 
LUN=1 

CALL  BTAKEC(LUN,  ISTAT)  ITAKE  CHARGE  OF  GPIB 

CALL  BDEVCL  ( IADDR ,  ISTAT)  ICLEAR  SELECTED  DEVICE 

CALL  ERRSET (64, ,  .  FALSE. ,  r  .  FALSE. , , )  IIGNORE  ERROR  MESSAGE  64 

CALL  BWRITE  (IADDR,  IADM,  4,  ISTAT)  DISPLAY  MODE 

AMM=10000. 

LENGTH1=35 
LENGTH2=21 
LENGTH3=80 
1053  FORM  AT  (F9. 5) 

AMN=0. 

WRITE  (5, 888) 

888  FORMAT G$DO  YOU  WANT  1)  IMPEDANCE  2) ADMITTANCE  ') 

READ  (5,  889)  ILOOP 

889  FORMAT  (12) 

WRITE  (5, 900) 

900  FORMAT (*$*, 'DO  YOU  WANT  1)LOW,  2) MED.  3) HIGH  SPEED  ') 

READ (5, 950)  IAVE 
950  FORMAT  (12) 


IPARALLEL  CIRCUIT  (ADMITTANCE) 

ISERIES  CIRCUIT  (IMPEDANCE) 

(TRIGGER  THE  READING 


20 


0103 

0104 

0105 

0106 

0107 

0106 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 
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10 

1000 


1060 

1070 

1111 

1020 

20 

1010 

1011 

1030 


1035 

1040 

1066 

1044 


WRITE  (5, 1000) 

FORMAT (/•$*, 'ENTER  THE  VOLTAGE  LEVEL  (IN  VOLTS):  ') 

READ(5f  1060,  ERR=10)  VOLTS 
ENCODE  <9,  1070,  A  (3))  VOLTS 
CALL  BWRITEOADDR,  A,  13,  ISTAT) 

FORMAT  (F10.0) 

FORMAT  (F9.4) 

NTIMES=1 

WRITE(5f  1111) 

FORMAT (  *5  ENTER  #  OF  POINTS  ') 

READ(5f  1020)  NPTS 
FORMAT  (14) 

NTIMES=NTIMES+1 

IF  (ILOOP  .EQ.  1)WRlTE(5f  1010) 

FORMAT  (/’$\ 'ENTER  FREQUENCY  (PARALLEL  RESONANCE)  (IN  KHZ):  ’) 

IF  (ILOOP  .EQ.  2)WRITE(5, 1011) 

FORMAT  (/’$*, ’ENTER  FREQUENCY  (SERIES  RESONANCE)  (IN  KHZ):  ’) 

READ  (5, 1060,  ERR=20)  START_FREQ 
ITIME=1 

WRITE  (5, 1030) 

FORMAT (’SENTER  STEP  FREQ  (IN  \ 

READ(5, 1060)STEP_FREQ 
IF(NTIMES  .LT.  5) THEN 
GOTO  3000 
ELSE 

TYPE  VDO  YOU  WANT  TO  SWEEP  FF 
READ  (5, 1035)  ISWEEP 
FORMAT  (A2) 

IF  (ISWEEP  .EQ.  'N*)THEN 
NTIMES=-10 
GO  TO  3000 
END  IF 
1=1 

ENCODE (8, 1066,  D(3)  )STEP_FREQ 
CALL  BWRlTEdADDR,  D,  12,  ISTAT) 

WRITE(5,  1040) 

FORMAT ( /  ,  '$ENTER  THE  STOPPING 
READ(5,  1060)  STOP 
ENCODE  (9, 1070,  C(3)) STOP 
FORMAT  (F8.4) 

CALL  BWRlTEdADDR,  C,  13,  ISTAT) 

ENCODE (9, 1070,  E(3)  )START_FREQ 
CALL  BWRlTEdADDR,  E,  13,  ISTAT) 

CALL  BWRlTEdADDR,  IAMM, 2,  ISTAT) 

CALL  BWRlTEdADDR,  IAEX,  2,  ISTAT) 

DO  J=1 , 80 
ANSWER  (J)=0 


KHZ)  :  *) 


=IEQ?’ 


FREQ  (KHZ)  ’) 


‘SPOT  FREQUENCY 


‘MANUAL  MODE 
ITRIGGER 
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0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 
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3000 


48 

49 


1050 


END  DO 

CALL  BREAD (IADDR,  ANSWER,  LENGTH3,  ISTAT) 

DECODE  ( 1 2, 1 050,  ANSWER  (5) )  RES  ( I ) 

DECODE  (12,1 050,  ANSWER  (21 ) )  XES  (I) 

DECODE  (9, 1 053,  ANSWER  (35) )  FREQ (I) 

IF  (FREQ(I)  .GE.  STOP)GOTO  88 

CALL  BWRITE  (IADDR,  IADV,  2,  ISTAT)  IADVANCE  FREQ. 

1=1+1 

GO  TO  1044 
END  IF 

STEPI=STEP_FREQ 

IF  (IAVE.EQ.1)  CALL  BWRITE  (IADDR,  IAV1 , 2,  ISTAT)  IAVERAGE  SPEED 
IF  GAVE  .EQ.  3)  CALL  BWRITE  (IADDR,  IAH1 , 2,  ISTAT)  IHIGH  SPEED 

IF  GAVE  .EQ.  2)  CALL  BWRITE  (IADDR,  IAHO,  2,  ISTAT)  INORMAL  SPEED 

CALL  BWRITE  (IADDR,  IAMM,  2,  ISTAT)  IMANUAL  MODE 

JFGLOOP  .EQ.  2)CALL  BWRITE  (IADDR,  IAC3, 2,  ISTAT)  ICIRCUIT  TYPE  3 
IFGLOOP  .EQ.  1)CALL  BWRITE  (IADDR,  IAC2, 2,  ISTAT)  ICIRCUIT  TYPE  2 

IF  (IT  .GT.  1)  CONTINUE 

DO  1=1,3 

IFGTIME  .GE.  2) THEN 

IF  (I  .EQ.  2)GOTO  85 

FREQUENCY=ABS(FR(I)  /( 1000.  *2*3. 14159) ) 

ENCODE  (9, 1070,  E  (3))  FREQUENCY  I  SPOT  STARTING  POINT 
CALL  BWRITE  ( IADDR,  E,  13,  ISTAT) 

FREQ(I)  =  FREQUENCY 

ELSE 

ENCODE  (9, 1070,  E  (3))  START_FREQ  I  SPOT  STARTING  POINT 
CALL  BWRITE  (IADDR,  E,  13,  ISTAT) 

FREQ ( I ) =START_FREQ 
END  IF 

CALL  BWRITE  (IADDR,  IAEX,  2,  ISTAT)  (TRIGGER  THE  READING 

DO  J=1 , 80 
ANSWER  (J)=0 
END  DO 
K=0 

CALL  BREADGADDR,  ANSWER,  LENGTH1 ,  ISTAT) 

DECODE  (12, 1050,  ANSWER  (5)  )  RES(I) 

DECODE  (12, 1050,  ANSWER  (21 )  )XES  (I) 

FORMAT  (El 2. 2) 

IFGTIME.  EQ.1.  AND. I.  EQ.  1)THEN 
FREQA=FREQ(1 ) 

XESA=XES (1 ) 

RESA=RES(1 ) 

ELSE 
END  IF 

FREQ(I)  =  FREQ(I)*1000 
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85  ARRAY  (1,1  )=RES(I) 

ARRAY  (l,2)=XES(l) 

ARRAY  (l,3)  =  1 . 

CONS  (l)=- (RES  (l)*RES(l)+XES(l)*XES(l) ) 

IF  (ITIME .  EQ .  1 )  START_FREQ=START_FREQ+STEP_FREQ 
END  DO 

C  SIMEQ  IS  A  SUBROUTINE  THAT  SOLVES  SIMELTANEOUS  EQUATIONS 

C  USED  TO  FIND  COEFICIENTS  OF  A  GENERALIZED  CIRCLE 

N=3  {NUMBER  OF  PARAMETERS 

CALL  SIMEQ(N,  ARRAY, CONS,  DEF) 

IF(N  .EQ.  0)GO  TO  20 

CENTERX=-DEF(1)/2.  !X  VALUE  OF  CENTER 

CENTER Y= -DEF (2)  /2 ,  !Y  VALUE  OF  CENTER 

RADIUS=  .  5*SQRT  ( DEF  ( 1 )  *DEF  ( 1 )  +DEF  (2)  *DEF  (2)  -4*DEF  (3) ) 

IF  (I  LOOP  .EQ.  2)  RESIS=  1 .  /  (2*RADIUS)  {RESISTANCE  (ADMITTANCE) 

IF  (I  LOOP  .EQ.  1 )  RESIS=2*RADIUS  {RESISTANCE  (IMPEDANCE) 

C  SLOPE  IS  THE  TAN  (THETA)  IN  EQUATION  1 

SLOPE1  =  ( (XES  ( 1 )  -CENTERY)  /  ( RES ( 1 ) -CENTERX+RADIUS) ) 

SLOPE2=  ( (XES (2)  -CENTERY)  /  (RES  (2) -CENTERX+RADIUS) ) 

SLOPE3=  ( (XES (3) -CENTERY)  /  (RES(3)-CENTERX+RADIUS) ) 

COEFF  (1,1 ) =SLOPE1  *RESIS*6 . 283185*FREQ  ( 1 ) 

COEFF  (2, 1  )=SLOPE2*RESIS*6.283l85*FREQ(2) 

COEFF  (1 , 2)  =  1 . 

COEFF(2,2)  =  1. 

COLUMN  ( 1 )  =  (2*3 . 1 41 59*FREQ  ( 1 ) )  **2 
COLUMN  (2)  =  (2*3. 14159*FREQ(2) )  **2 

N=2 

CALL  SIMEQ (N,  COEFF,  COLUMN,  ANS) 

IF (N  .EQ.  0) GOTO  20 

C  ANS(1 )  is  the  Inductor 

C  ANS(2)  is  the  Resonant  Freq. 

IF  (ANS  (2)  .LT.  0)GOTO  20 

ANS  ( 1 )  =  1  /ANS  ( 1 ) 

ANS(2)  =  SQRT (ANS  (2) ) 

IF  (ANS  (2)  .EQ.  0)  GOTO  20 
RESONANCE=ANS(2)  /6283. 185 
WRITE  (5,  77)  ANS (2)  /6283 . 1 85 
77  FORMAT (*  FREQ=  \F10.3) 

IF  (ABS  ( ANS  (2)  -OMEGA)  /6283 ,  LE .  ( .  0005*FREQA)  .OR.  ITIME. GE.  25)  THEN 
ANS  (2)  =  (OMEG A+ANS  (2) )  /2 . 

ANS  ( 1 )  =  ( Rl N DUCTOR+ANS  ( 1 ) )  /2 . 

GOTO  50 
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ELSE 

0MEGA=ANS(2) 

RINDUCTOR=ANS  (1 ) 

END  IF 

C  ANGLES  OF  OTHER  2  POINTS 

TRANSX1=RES(2)-CENTERX 
TRANSY1=XES  (2)  -CENTERY 
TRANSA1  =ATAN2 (TRANSY 1  f  TRANSX1 ) 

X2=COS  (TRANSA1+2 . 094395)  *RADIUS+CENTERX 
X3=COS  (TRANSA1-2 . 094395)  *RADIUS+CENTERX 
Y2=SIN  (TRANSA1  +2 . 094395)  *RADIUS+CENTERY 
Y3=SIN  (TRANSA1  -2 . 094395)  *RADIUS+CENTERY 

SLOPE  1=  (Y2-CENTERY)  /  (X2-CENTERX+RADIUS) 

SLOPE3=  (Y3-CENTERY)  /  (X3-CENTERX+RADIUS) 
B1=-SLOPE1*RESIS 
B2=-SLOPE2*RESIS 
B3=-SLOPE3*RESIS 

FR  ( 1 )  =-B1-SQRT(B1  *B1+4 .  *ANS  ( 1 )  **2*  (ANS  (2)*ANS  (2)  ) ) 

FR  ( 1 )  =FR  ( 1 )  /  (2*ANS  ( 1 ) ) 

FR  (2)  =FREQ  (2)  *2*3 .14159 

FR  (3 )  =“B3“SQRT ( B3*B3+4 .  *ANS  ( 1 ) **2*  (ANS (2)  *ANS  (2) ) ) 
FR(3)=FR(3)/(2*ANS(D) 

ITIME=ITIME+1 

IFCITIME  .GE.  2) GOTO  49 
50  R=SQRT  (RESA*RESA+XESA*XESA) 

THETA=ATAN2 (XESA,  RESA) 

AO=SQRT  ( (CENTERX-RESA)  **2  +  (CENTER  Y-XESA)  **2) 

CON  VERT=  1 . 74532E-02 
TIM=359.  /NPTS 
Pl=3. 14159 

TRANSX=RESA-CENTERX 
TRANSY=XESA-CENTERY 
THETA=ATAN2  (TRANSY,  TRANSX) 

1=1 

DO  IT=1 ,  NPTS 

DELTA2=IT*TIM*C0NVERT+THETA 
X3=  (COS(DELTA2)  *RADIUS)+CENTERX 
Y3=  (SIN  (DELTA2)  *RADIUS)  +CENTERY 

SLOPE3=  (Y3-CENTERY)  /  (X3-CENTERX+RADIUS) 

B3=-SLOPE3*RESIS 

FR  E  ( I )  =-B3+SQRT  ( B3*B3+4 .  *ANS  ( 1 )  **2*  (ANS  (2 )  *ANS  (2 ) ) ) 
FRE(I)  =  FRE(I)  / (2.  *ANS (1) ) 

FR EQU ENCY= ABS  ( FR E  ( I )  /  ( 1 000 .  *2*3 . 1 41 59) ) 

9  TOLERANCE1=2 . 5*RESONANCE 

IF  (  FREQUENCY  ,GT,  TOLERANCE  1) GOTO  70 


24 


0307 

0306 

0306 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

0340 

0341 

0342 

0343 

0344 

0345 

0346 

0347 

0348 

0349 

0350 

0351 

0352 

0353 

0354 

0355 

0356 

0357 


NRL  MEMO  REPORT  5384 


ENCODE  (9f  1070,  E  (3))  FREQUENCY  I  SPOT  STARTING  POINT 
CALL  BWRITEOADDR,  E,  13,  ISTAT) 


CALL  BWRITEOADDR,  IAEX,  2,  ISTAT) 


FTRIGGER  THE  READING 


DO  J=1 , 80 

ANSWER  (J)=0 
END  DO 
K=0 

CALL  BREAD (IADDR,  ANSWER,  LENGTH  1,  ISTAT) 

DECODE ( 12, 1050,  ANSWER  (5) )  RES(I) 

DECODE  (12,1050,  ANS WE R  ( 2 1 ) )  XES  ( I ) 

IF(RES(I)  .GT.  1.3E6  .OR.  XES(IT)  .GT.  1.3E6)GOT0  70 
FR  EQ ( I ) = FR  EQU ENCY 

MAG  (l)=SQRT(RES(l)*RES(l)+XES(l)*XES(l) ) 

1=1+1 
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END  DO 

CALL  RSORT (RES,  FREQ,  XES,  1-1, 1) 


GS=RES  (1 ) 

BS=XES ( 1 ) 

FS=FREQ(1) 

TAN0=ATAN2(BS,GS) 

CALL  RSORT(FREQ,  RES,  XES,  1-1,1) 

FLAG=0 

IT=2 

XESMAX=-1E12 

RESMAX=-1E12 

XESMIN=1E12 

RESMIN=1E12 

ADMMIN=1E12 

ADMMAX=0 

TANMIN=1E12 

DO  J=2, 1-1 

IF (XES(J)  .GT.  BS  .AND.  FLAG  .EQ.O  .AND.  XES(J)  .GT.  0)THEN 
FREQ(1)=FREQ(J) 

RES(1)=RES(J) 

XES(  1)=XES(J) 

GO  TO  500 
ELSE 

FREQ(IT)  =  FREQ(J) 

RES  (IT)=RES  (J) 

XES  (IT)=XES  (J) 

TAN1=ATAN2(XES  (IT) ,  RES  (IT) ) 

FLAG=1 

ADM(IT)=SQRT(RES(IT)**2  +XES(IT)**2) 

IF(FREQOT)  .LE.  FS  .AND.  IT  .GT.  1)THEN 
IF  (XES(IT)  .GE.  XESMAX)  THEN 
XESMAX=XES  (IT) 

FH=FREQ(IT) 

GH=RES(IT) 

BH=XES  (IT) 

END  IF 

IF(ADM(IT)  .GT.  ADMMAX  .AND.  ILOOP  .EQ.  2)THEN 
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500 

1109 

1110 

1112 


ADMMAX=ADM  (IT) 

FM=FREQ(IT) 

GM=RES(IT) 

BM=XES(IT) 

END  IF 

IF(XESUT-I)  .LT.  0  .AND.  XES(IT)  .GT.  0)THEN 
FRES=FREQ(IT) 

GR=RES(IT) 

BR=XES(IT) 

END  IF 
END  IF 

IF(FREQGT)  .GE.  FS  .AND.  IT  .GT.  1)THEN 
IF (XES(IT)  .LE.  XESMIN)THEN 
XESMIN=XES  (IT) 

FL=FREQ(IT) 

GL=RES(IT) 

BL=XES  (IT) 

END  IF 

IF (XESOT-1 )  .LT.  0  .AND.  XES(IT)  .GT.  0)THEN 
FRES=FREQ(IT) 

GR=RES(IT) 

BR=XES(IT) 

END  IF 

IF  (ADM  (IT)  .GT.  ADMMAX  .AND.  ILOOP  .EQ.  1)THEN 
ADMMAX=ADM  (IT) 

FN=FREQ(IT) 

GN=RES(IT) 

BN=XES(IT) 

ENDIF 
END  IF 
IT=IT+1 
END  IF 
END  DO 

IF  ( ILOOP.  EQ.1)  THEN 
WRITE  (5, 1109) 

FORMAT(///,  9X,TYPE',11X,  'FREQUENCY  (KHZ)  \  2X, 'RESISTANCE',  2X, 

1 'REACTANCE') 

ELSE 

WRITE(5f  1110) 

FORMAT  ( /  /  / ,  9X,  'TYPE  \1  IX,  'FREQUENCY  (KHZ)',  2X,  'CONDUCTANCE',  2X, 
1'SUSCEPTANCE') 

END  IF 

WRITE  (5, 1112) 

FORMAT (9X, ' - ',  11X, ' - \2X,' - V 

12X,  * - ',/) 

IF(FRES  .EQ.  0)THEN 
IF(BL  .LT.  0)THEN 
FRES=FS 
GR=GS 
BR=BS 

WRITE(5, 1130)FRES,  GR,  BR 
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1120 

1130 

1131 


1270 


ELSE 

WRITE(5, 1120) 

FORMAT  ( IX,  "RESONANCE',  1 7Xt '  DOES  NOT  EXIST') 

END  IF 
ELSE 

IF  (I  LOOP  .EQ.  2)WRITE(5,  1 130)FRES,  GR,  BR 

FORMAT (IX, 'RESONANCE',  16X,F8.3,7X,  E8.3, 5X,  E9. 3) 

IFOLOOP  .EQ.  1 ) WRITE (5, 1 131  )FRES,  GR,  BR 

FORMAT (IX,  'ANTI RESONANCE',  12X,  F8.3,  7X,  E8. 3, 5X,  E9.3) 

END  IF 

CALL  BWRITEOADDR,  IAMM,  2,  ISTAT)  IMANUAL  MODE 

FLAG1=0 

ENCODE  (8, 1066,  D  (3)  )  .  1  ISTEP  SIZE 

CALL  BWRITEOADDR,  D,  12,  ISTAT) 

IFOLOOP  .EQ.  1)CALL  BWRITEOADDR,  IAC3, 2,  ISTAT) 

IFOLOOP  .EQ.  2)  CALL  BWRITEOADDR,  IAC2, 2,  ISTAT) 

IFOLOOP  .EQ.  2) ENCODE (9,  1070,  E(3)) FL 
IF  (I LOOP  .  EQ .  1 )  ENCODE (9, 1070,  E<3) ) FH 

CALL  BWRITEOADDR,  E,  13,  ISTAT)  ISPOT 

GP=0 

CALL  BWRITEOADDR,  I AEX, 2,  ISTAT) 

DO  J=1 , 80 
ANSWER  (J)=0 
END  DO 

CALL  BREAD (IADDR,  ANSWER,  LENGTHS,  ISTAT) 

DECODE(  12, 1050,  ANSWER  (5) )  A1 
DECODE  ( 12, 1050,  ANSWER  (21 ) )  B1 
DECODE  (9, 1053,  ANSWER  (35) )  FI 
IF (A1  .GE.  GP  )THEN 
FP=F1 
GP=A1 
BP=B1 

IFOLOOP  .EQ.  2) CALL  BWRITEOADDR,  IADV, 2,  ISTAT)  (ADVANCE 

IFOLOOP  .EQ.  1)CALL  BWRITEOADDR,  IREV, 2,  ISTAT)  (REVERSE 

GO  TO  1270 
ELSE 

FLAG1=FLAG1+1 
IF(FLAG1  .EQ.  1)  THEN 

IF  (ILOOP .  EQ .  2) CALL  BWRITE  (IADDR,  IADV,  2,  ISTAT) 

IFOLOOP.  EQ.1) CALL  BWRITEOADDR,  IREV,  2,  ISTAT) 

FP=F1 
GP=A1 
BP=B1 
GOTO  1270 
END  IF 
END  IF 
NUM=8 
CRX(1)=GS 
CRY(1  )=BS 
CRX(2)=GP 
CRY  (2)  =  BP 
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CRX(3)=QM 
CRY  (3)=BM 
CRX(4)=GN 
CRY  (4)=BN 
CRX(5)=GH 
CRY  (5)=BH 
CRX(6)=GL 
CRY  (6)=BL 

IF(FRES  .EQ.  0)THEN 
NUM=NUM-1 
ELSE 

CRX(7)=GR 
CRY  (7)=BR 
END  IF 

IF  (FA  .EQ.  0)THEN 
NUM=NUM-1 
ELSE 

CRX(8)=GA 
CRY  (8)=BA 
END  IF 

IFGLOOP  .EQ.  2)WRITE(5, 1 160)FS,  GS,  BS 
IF (ILOOP  .EQ.  1 ) WRITE (5f  1 160) FP, GP,  BP 

FORMAT (IX,  ’SERIES  RESONANCE’,  9X,  F8.3,  7X,  E8.3,5X,  E9.3) 
IFGLOOP  .EQ.  2) WRITE (5, 1170) FP, QP,  BP 
IFGLOOP  .EQ.  1)WRrTE(5, 1 170)FS, GS,  BS 

FORMAT (IX,  ’PARALLEL  RESONANCE’,  7X,  F8.3,  7X,  E8.3.5X,  E9.3) 
IFGLOOP  .EQ.  2)WRITE(5,  1 180)FM,  GM,  BM 

FORMAT (IX,  'MAXIMUM  ADMITTANCE’,  7X,  F8.3,  7X,  E8.3, 5X,  E9.3) 
IFGLOOP  .EQ.  1)  WRITE (5, 1280) FN, GN,  BN 

FORMAT (IX,  'MAXIMUM  IMPEDANCE’,  8X,  F8. 3,  7X,  E8.3, 5X,  E9.3) 
WRITE  (5,  1190)  FH,GH,BH 

FORMAT (IX,  ’HIGH  POINT,  15X,  F8.3,  7X,  E8. 3, 5X,  E9.3) 

WRITE (5, 1200)  FL,  GL,  BL 

FORMAT (IX,  'LOW  POINT,  16X,  F8.3,  7X,  E8.3,5X,  E9.3) 

WRITE  (5, 1300)7 

FORMAT(//,’  ’,  A1 ) 

IFGLOOP  .EQ.  2) TYPE  *, 'ADMITTANCE  LOOP* 

IFGLOOP  .EQ.  1)TYPE  *,  'IMPEDANCE  LOOP 
CALL  FILE  (IT-1,  RES,  XES) 

TYPE  *,'  ’ 

TYPE  *,  ’CRITICAL  FREQUENCIES’ 

CALL  FILE  (NUM,CRX,  CRY) 

TYPE  V  * 

TYPE  *, ’FREQ  VS.  MAG.' 

CALL  FILEGT-1,  FREQ,  MAG) 

CALL  BDEVCL  ( IADDR ,  ISTAT )  ICLEAR  SELECTED 

END 
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Running  MI  CAM 


RUN  SUPER 

DO  YOU  WANT  1 )  IMPEDANCE  2  )  ADMITTANCE  1 
DO  YOU  WANT  1)  LOW,  2)  MED,  3)  HIGH  SPEED  3 

ENTER  THE  VOLTAGE  LEVEL  (  IN  VOLTS  ):  .1 

ENTER  #  OF  POINTS  100 

ENTER  FREQUENCY  ( PARALLEL  RESONANCE )  (  IN  KHZ ) :  15 

ENTER  STEP  FREQ  (  IN  KHZ  ) :  .5 

FREQ-  15.304 

FREQ-  15.293 

FREQ-  15.291 

TYPE 


ANTIRESONANCE 
SERIES  RESONANCE 
PARALLEL  RESONANCE 
MAXIMUM  IMPEDANCE 
HIGH  POINT 
LOW  POINT 

IMPEDANCE  LOOP 
FILE  NAME  FOR  X,Y  DATA:  DATA1 

CRITICAL  FREQUENCIES 

FILE  NAME  FOR  X,Y  DATA:  DATA2 

FREQ  VS.  MAG. 

FILE  NAME  FOR  X,Y  DATA:  DATA3 

$ 


FREQUENCY (  KHZ )  RESISTANCE 


15.291 

. 107E+O7 

14.384 

. 485E-02 

15.291 

. 107E+07 

15.291 

.107E+07 

15.284 

.530E+06 

15.297 

.630E+06 
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REACTANCE 

0 . 100E405 
-.288E-02 
0 . 100E+05 
0 . 100E+05 
0 . 530E+06 
-.540E+06 
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RUN  SUPER 

DO  YOU  WANT  1 )  IMPEDANCE  2  )  ADMITTANCE  2 
DO  YOU  WANT  1)  LOW,  2)  MED,  3)  HIGH  SPEED  3 

ENTER  THE  VOLTAGE  LEVEL  (IN  VOLTS):  .1 

ENTER  #  OF  POINTS  100 

ENTER  FREQUENCY (  SERIES  RESONANCE)  (IN  KZH ) :  14 

ENTER  STEP  FREQ  (IN  KHZ):  .5 

FREQ*  14.390 
FREQ*  14.372 
FREQ*  14.372 

TYPE 


RESONANCE 
SERIES  RESONANCE 
PARALLEL  RESONANCE 
MAXIMUM  ADMITTANCE 
HIGH  POINT 
LOW  POINT 

ADMITTANCE  LOOP 
FILE  NAME  FOR  X,Y  DATA;  DATA1 

CRITICAL  FREQUENCIES 

FILE  NAME  FOR  X,Y  DATA;  DATA2 

FREQ  VS.  MAG. 

FILE  NAME  FOR  X,Y  DATA:  DATA3 

$ 


FREQUENCY ( KHZ ) 

CONDUCTANCE 

SUSCEPTANCE 

14.375 

.678E-02 

0 . 300E-04 

14.373 

.679E-02 

0 . 540E—03 

15.291 

.  107E4-07 

0 . OOOE+OO 

14.373 

.677E-02 

0 . 790E—03 

14.357 

.311E-02 

0 . 379E—02 

14.391 

.327E-02 

— . 316E—02 
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